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INTRODUCTION 


This  presentation  will  describe  our  initial  experience  with  the 
measurement  and  interpretation  of  the  two-wavevector-f requency 
spectra  of  nonhomogeneous  but  stationary  fields.  Let  us  begin  with 
the  definitions  appropriate  to  the  nonhomogeneous  stationary  case. 


SLIDE  1 

As  shown  at  the  top  of  this  slide,  a  space-varying  spectrum  is 
defined  as  the  Fourier  transform  of  the  correlation  of  a 
nonhomogeneous,  stationary  field  variable  that  is  a  function  of  an 
absolute  spatial  location,  x,  of  one  observation  of  the  field 
variable,  rhe  spatial  separation,  ,  between  observations,  and  the 
time  dela'',  T,  betv/een  observations. 

Recall,  from  Dr.  Strawderman ' s  earlier  talk,  that  the  two-wavevector- 
frequency  spectrum  of  a  nonhomogeneous,  stationary  field  variable  is 
obtained  by  taking  a  Fourier  transform  of  the  space-varying  spectrum, 
as  shown  at  the  bottom  of  the  slide.  The  wavevector,  n,  is  the 
Fourier  conjugate  variable  of  the  absolute  spatial  vector,  x;  the 
wavevector,  k,  is  the  conjugate  variable  of  the  spatial  separation 
vector,  f;  and  the  frequency,  o,  is  the  conjugate  variable  of  the 
time  delay,  T.  Note  that  the  number  of  the  wavevectors  and  their 
dimensions  must  exactly  match  the  number  of  spatial  variables  and 
their  dimensions  that  are  present  in  the  correlation  field.  Spatial 
spectra  resulting  from  nonhomogeneous  fields  must  necessarily  have  a 
higher  dimensionality  than  its  homogeneous  counterpart  because  the 
correlation  field  is  a  function  of  both  a  location  variable  and  a 
separation  variable. 
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PROJECTOR  OFF 

Experimentally  this  increase  in  the  dimensionality  of  the  problem 
represents  a  significant  complication,  for  two  reasons.  First,  the 
process  of  viewing  wavevector-f requency  spectra  is  difficult  to  begin 
with  because  the  surface,  in  general,  exceeds  three  dimensions  for 
homogeneous  fields.  In  the  case  of  a  stationary  nonhomogeneous  field 
one  is  faced  with  the  task  of  examining  an  eight-dimensional  figure 
unless  the  physics  of  the  problem  reduce  the  number  of  nonhomogeneous 
directions.  The  second  difficulty  is  the  lack  of  experience  with  the 
wavevector  spectra  associated  with  the  length  variable  characterizing 
the  nonhomogeneous  nature  of  the  problem.  Thus,  even  after  coming  to 
grips  with  the  analytical  and  numerical  problems  of  the  spectral 
estimation,  we  are  faced  with  such  data  analysis  problems  as  "what 
does  it  look  like?"  and  "is  it  useful?"  or  more  exactly  "was  it  worth 
doing?"  In  this  presentation  we  examine  these  questions  for  an 
inherently  low  dimensional  structure;  that  is,  a  linear  hydrophone 
array  in  near  axial  flow. 


SLIDE  2 

The  previous  slide  illustrated  the  conceptual  form  of  the  desired 
quantity.  However,  other  forms  are  available  and  are  important  when 
developing  spectral  estimation  algorithms,  as  will  be  shown.  In  this 
slide  we  show  a  computational  form  as  well  as  the  preferred 
conceptual  form.  The  computational  form  is  seen  to  be  nested  Fourier 
transforms  of  the  fundamental  spatial  variables  of  the  problem.  The 
second  equation  gives  the  computational  form  in  terms  of  the 
experimentally  important  quantity;  cross-spectral  density.  However, 
the  resulting  wavevector  does  not  maintain  its  familiar  relation  to 
the  correlation  separation  variable  and  thus  eliminates  decades  of 
experience.  The  preferred  conceptual  form  maintains  this  relation; 
thus,  we  use  the  computational  form  initially  and  then  by  a 
transformation  complete  the  calculation  with  the  conceptual  form. 
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Forms  of  Two-Wavevector-Frequency  Spectra 
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SLIDE  3 

In  much  of  our  experimental  work  our  signal-processing  equipment  is 
designed  to  efficiently  calculate  the  matrix  of  cross-spectral 
densities  associated  with  an  experiment  containing  n  transducers. 

Thus  the  cross-spectral  density  (CSD)  matrix  is  a  good  starting  point 
for  a  tv;o-wavevector  calculation.  This  slide  shows  the  computational 
form  of  the  two-wavevector- frequency  spectra  in  terms  of  the  CSD 
matrix  components  for  a  one-dimensional  structure. 
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Calculation  of  Two-Wavenumber  Spectra 
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SLIDE  4 

If  we  solve  the  previous  equation  at  the  unique  set  of  v;avenumbers 
associated  with  the  fast  Fourier  transfer  algorithm,  then  the 
resulting  two-wavevector-f requency  spectral  estimation  algorithm  is 
seen  to  be  a  pair  of  nested  FFTs.  A  final  transformation  yields  the 
preferred  form. 
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Evaluate  Ka  at  the  FFT  Wavenumbers 
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SLIDE  5 

In  the  remainder  of  this  presentation  we  will  show  examples  of  one- 
and  two-wavenumber-frequency  spectra.  First,  canonical  forms  of 
correlation  fields  will  be  used  to  illustrate  fundamental 
"footprints"  that  will  also  be  observed  in  the  measured  data.  This 
step-by-step  approach  v/as  necessary  because  of  lack  of  experience 
v/ith  the  subject  density  function.  In  this  slide  we  show  the 
simplest  field  possible;  that  is,  a  homogeneous  one  for  which  the  CSD 
matrix  is  just  the  identity  matrix.  The  k  wavenumber  is  associated 
with  the  separation  variable,  and  the  wavenumber  is  associated  v/ith 
the  location  variable.  Because  the  field  is  homogeneous,  there  is  no 
energy  in  the  dimension;  and  because  the  field  is  isotropic,  the 
energy  in  the  k  dimension  is  wavenumber  white.  Note  that  at  all 
nonzero  values  of  ^  the  spectral  values  result  from  "side  lobe 
leakage"  in  the  estimation  algorithm. 
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SLIDE  6 

In  this  slide  we  show  a  fundamental  pattern.  The  figure  is  a  contour 
plot  of  the  three-dimensional  space  associated  with  the  calculated 
two-wavenumber- frequency  spectra.  The  contours  are  constant  values 
of  a  pressure  density  function  with  the  dimensions  /xPa^Hz  (rad/m)  2. 

The  spectrum  shown  here  results  from  a  unidirectional  traveling  wave 
decaying  exponentially  in  space.  The  peak  in  the  k  dimension  is  the 
familiar  propagation  wavenumber.  Because  the  field  is  now 
nonhomogeneous ,  there  is  a  distinctive  distribution  of  energy  in  the 
^  dimension. 
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SLIDE  7 

Before  showing  measured  two-waveni,  spectra,  it  is  useful  to 

examine  the  more  familiar  one-wavenumber-frequency  spectral  densiry 
of  the  m.easurement .  This  is  equivalent  to  the  spaced  average  spectra 
discussed  by  Dr.  Strawderman.  Recall  that  the  structure  is  a  linear 
array  of  equally  spaced  hydrophones  in  axial  flow.  The  numerical 
derails  of  the  experiment  are  not  important  to  this  presentation,  nor 
does  time  allow  a  detailed  description.  To  orient  you  to  the  figure, 
note  the  obvious  "sonic"  cone  that  separates  the  acoustic  from 
ncnacoustic  v/avenumbers .  Significant  amounts  of  energy  are  present 
in  both  domains.  Ilext  please  note  the  narrow  band  of  energy  that 
extends  across  all  wavenumbers  at  9.5  Hz.  This  energy  results  from  a 
supersonic  longitudinal  vibration  in  the  structure.  The  structure 
consists  of  a  c.hain  of  long  sections  that  creates  a  periodic 
mechanical  discontinuity  at  the  juncture  of  the  sections.  The 
measured  pressure  field  contains  "hot  spots"  at  each  of  the 
junctures,  creating  a  distinctly  nonhomogeneous  pressure  field.  T.his 
periodic  .modulation  of  the  amplitude  of  the  pressure  field  manifests 
itself  in  the  one-wavenumber  spectra  as  periodic  "sidebands”  of 
energy  •,;ith  a  separation  proportional  to  the  reciprocal  of  the 
lengths  of  the  sections.  The  wavenumber  structure  is  quite 
complicated,  and  only  rather  obvious  qualitative  infor.mation  can  be 
realized  from  this  figure. 


SLIDE  8 


This  slide  illustrates  the  two-wavenumber-frequency  spectra  of  the 
measurement  at  a  single  frequency;  that  is,  9.5  Hz.  The  horizontal 
axis  is  the  homogeneous  wavenumber  k,  and  the  vertical  axis  is  the 
nonhomogeneous  wavenumber  n.  One  immediately  sees  the  "footprint"  of 
an  exponentially  decaying  longitudinal  wave.  Quantitative  analysis 
of  the  ^  dependence  of  the  energy  is  a  measure  of  the  damping  factor 
of  the  wave.  The  velocity  of  the  wave  is  obtained  from  the  peak  in 
the  pattern  at  m  equal  to  zero.  Note  also  the  apparent  lack  of  a 
homogeneous  component  to  the  field  as  evident  from  the  lack  of  a 
ridge  in  k  at  zero  n. 
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SLIDE  9 

In  this  slide  we  show  the  two-wavenumber  spectra  of  the  data  at  a 
frequency  of  79  Hz.  This  frequency  was  chosen  because  it  contained 
little  energy  from  the  previously  examined  vibrational  v/ave  and  also 
because  it  contained  significant  sonic  energy.  Here  we  see  clear 
evidence  of  homogeneous  acoustic  energy  and  some  homogeneous 
nonacoustic  waves.  However,  the  general  nonacoustic  energy  is, 
strangely,  uniformly  nonhomogeneous. 


18 


TD  8429 


SLIDE  10 

In  summary,  we  have  shown  examples  of  two-wavenumber  spectra 
calculated  from  both  theoretical  and  measured  cross-spectral  density 
matrices.  Canonical  forms  of  common  fields  were  used  to  establish 
identifiable  patterns  that  were  then  observed  in  the  data.  In 
particular,  the  approach  appears  to  be  useful  in  the  quantitative 
identification  of  nonhomogeneous ,  exponentially  damped  waves. 
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